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Abstract 

“Arithmetic random waves” are the Gaussian Laplace eigenfunctions on the two- 
dimensional torus [RWl IKKW] . In this paper we find that their nodal length converges 
to a non-universal (non-Gaussian) limiting distribution, depending on the angular dis¬ 
tribution of lattice points lying on circles. 

Our argument has two main ingredients. An explicit derivation of the Wiener-Ito 
chaos expansion for the nodal length shows that it is dominated by its 4th order chaos 
component (in particular, somewhat surprisingly, the second order chaos component 
vanishes). The rest of the argument relies on the precise analysis of the fourth order 
chaotic component. 

Keywords and Phrases: Arithmetic Random Waves, Nodal Lines, Non-Central Limit 
Theorem, Berry’s Cancellation. 

AMS Classification: 60G60, 60D05, 60B10, 58J50, 35P20 


1 Introduction and main results 


1.1 Arithmetic random waves 


Let T := 'S? jT? be the standard 2-torus and A the Laplacian on T. We are interested in 
the (totally discrete) spectrum of A i.e., eigenvalues Li > 0 of the Schrodinger equation 

A/ + i?/ = 0. (1.1) 

Let 

S = {n ■. n = + for some a, 6 G Z} 

be the collection of all numbers expressible as a sum of two squares. Then, the eigenvalues 
of (11.11) (also called energy levels of the torus) are all numbers of the form En = Ivr^n with 
n G S. 

In order to describe the Laplace eigenspace corresponding to En, denote by A„ the set of 
frequencies: 

An := {A = (Ai, A2) G Z^ : A^ -|- A2 = n} 

whose cardinality 

Vn := |A„| = r2(n) (1.2) 

equals the number of ways to express n as a sum of two squares. (Geometrically, An is the 
collection of all standard lattice points lying on the centred circle with radius y/n.) For 
X G An denote the complex exponential associated to the frequency A 


e\(x) = exp(27ri(A, x)) 
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with X = {xi,X 2 ) € T. The collection 


{eA(x)}AeA„ 


of the complex exponentials corresponding to the frequencies A G A„, is an L^-orthonormal 
basis of the eigenspace £n of A corresponding to the eigenvalue En- In particular, the 
dimension of £n is 

dimT^j — A/^ — |A-n| 

(cf. (I 1 . 2 p i. The number Mn is subject to large and erratic fluctuations; it grows |La] on 
average as Vlog n, but could be as small as 8 for (an infinite sequence of) prime numbers 
p = 1 mod 4, or as large as a power of log n. 

Following [RW| and |KKW] , we define the arithmetic random waves (also called random 
Gaussian toral Laplace eigenfunctions) to be the random fields 


Tn{x) 



a\ex{x), 

AG Att, 


X G T, 


(1.3) 


where the coefficients ax are standard complex-Gaussian random variables verifying the 
following properties: ax is stochastically independent of whenever 7 ^ {A, —A}, and 


O-A = ax 


(ensuring that the T„ are real-valued)Ji] By the definition (jl.3p . Tn is a stationary (i.e. the 
law of Tn is invariant under all the translations 

/(•) fix' + ■), 

x' G T), centered Gaussian random field with covariance function 

r„(x,x') = rnix-x) := E[T„(x)r„(x')] = ^ exix-x) = cos (27r(x - x', A)) , 


AgA^j 


AGAtj 


x,x' G T (by the standard abuse of notation for stationary fields). Note that r„(0) = 1, i.e. 
Tn has unit variance. 


1.2 Nodal length: mean and variance 

Consider the total nodal length of the random eigenfunctions, i.e. the collection {Cn}n&s of 
all random variables with the form 

Cn := length(T-Ho}). (1.4) 

The expected value of Cn was computed in |RWj to be 

E[T„] = (1-5) 

consistent with Yau’s conjecture [YajlDF] . The more challenging question of the asymptotic 
behaviour of the variance Var(£„) of Cn was addressed in [RW] . and fully resolved in |KKW] 
as follows. 

^From now on, we assume that every random object considered in this paper is defined on a common 
probability space (tl, Y, P), with E denoting mathematical expectation with respect to P. 
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Given n £ S, define a probability measure /i^ on the unit circle 5^ C supported on 
angles corresponding to lattice points in A„: 

1 

IJ'n ■ = 


AA 


V 6^. 

AeAn 


It is known m that for a density 1 sequence of numbers {rij} C S the angles of lattice 
points in A„ tend to be equidistributed, in the sense that 

d(f> 


^ 27r 


( 1 . 6 ) 


(where ^ indicates weak-* convergence of probability measures, and dcf) stands for the 
Lebesgue measure on 5^). However the sequence {fJ,n}nes has other weak-* adherent points 
[CH IKKW] (called attainable measures ), partially classified in |KWj . 

It was proved in |KKW] that one has 


where 


Cn. — 


and, for a measure ^ on 5^, 

Kk) = 

are the Fourier coefficients of /U on the unit circle. As 

|/J'n(4)| < 1 


^■^(1 + OAA^oo(l)), 

(1.7) 

1 + /^n(4)^ 

512 

(1.8) 

z~^ dpL{z), k £ Z, 



by the triangle inequality, the result (|1.7p shows that the true order of magnitude of Var(/1„) 
is this is of smaller order than what would be a natural guess, namely this situation 
(customarily called arithmetic Berry’s cancellation, see [KKW] 1 is similar to the cancellation 
phenomenon observed by Berry in a different setting, see [EillM]. 

In addition, dni) shows that, in order for Var(£„) to exhibit an asymptotic law (equivalent 
to {cn} in (II.8p being convergent along a subsequence) we need to pass to a subsequence 
{uj} C S such that the limit 

lim I/In, (4) I 

exists. For example, if {%■} C 5 is a subsequence such that /^n, h for some probability 
measure /x on 5^, then (II.7p reads (under the usual extra-assumption oo) 

Var(£n,) ~ (1-9) 


with 


c{ia) 


512 


where, here and for the rest of the paper, we write On ~ bn to indicate that the two 
positive sequences {on} and {bn} are such that anjbn —)• 1, as n —)• oo. Here, the set of 
the possible values for the 4th Fourier coefficient /I(4) attains the whole interval [—1,1] 
(see |KKW( IKWj l. This implies in particular that the possible values of the asymptotic 
constant c(//) attain the whole interval [^, ; the above is a complete classification of 

the asymptotic behaviour of Var(/ln). 
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1.3 Statement of the main results: asymptotic distribution of the nodal 
length 


Our main goal is the study of the fine asymptotic behaviour, as Nn —>■ oo, of the distributions 
of the sequence of normalised random variables 


Cn - ^[Cn] 
\/Var(£„) ’ 


n G S, 


( 1 . 10 ) 


(this is equivalent to studying Cn^ along subsequences C S satisfying ^ oo; 

note that it is possible to choose a full density subsequence in S as above). Since the variance 
(|1.7h diverges to inhnity, it seems reasonable to expect a central limit result, that is, that 
the sequence {Cn} converges in distribution to a standard Gaussian random variable. Our 
findings not only contradict this (somewhat naive) prediction, but also classify all the weak- 
* adherent points of the probability distributions associated with the collection of random 
variables : n G S'! (where the adherent points are in the sense of weak-* convergence 
of probability measures). In particular, we will show that such a set of weak-* adherent 
points coincides with the collection of probability distributions associated with a family 
of linear combinations of two independent squared Gaussian random variables; these linear 
combinations are parameterized by the adherent points of the sequence {|/4i(4)|} of real non¬ 
negative numbers < 1. This will show the remarkable fact that the angular distribution of 
A„ (or, more specihcally, the 4th Fourier coefficient of fin) does not only prescribe the leading 
term of the nodal length variance Var(£„), but, in addition, it prescribes the asymptotic 
distribution of Cn- 


To state our results formally, we will need some more notation. For r] G [0,1], let A4,y be 
the random variable 

■= 

2vl -F 

where X = {Xi,X 2 ) are independent standard Gaussians. Note that for rji ^ rj 2 the 
distributions of Mrn and are genuinely different; this follows for example from the 
observation that the support of the distribution of M-rj is 



Our first main result establishes a limiting law for the nodal length distribution for 
subsequences {nj}j>i C S provided that the numerical sequence 


{|/In, (4)| : j > 1} 


of non-negative numbers is convergent. As it was mentioned above, for some full density 
subsequence {nj}j>i C S the corresponding lattice points A„^. are asymptotically equidis- 
tributed (jl.6h . so that for this subsequence, in particular. 


7‘nj(4) —)■ 0. 

More generally, if for some subsequence {nj}j>i C S the angular distribution of the corre¬ 
sponding lattice points converges to fi, i.e. fin^ fi, where fi is some probability measure 
on 5^, then 

/ 2 „.( 4 ) ^ / 2 ( 4 ). 
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From now on, we use the symbol —?• to denote convergence in distribution of random 
variables; similarly, we will write X = y to indicate that the random variables X and Y have 
the same distribution. Observe that a sequence of random variables converges in distribution 
if and only if the corresponding sequence of probability laws is weak-* convergent. We shall 
however use the sentence “convergence in distribution” (resp. “weak-* convergence”) for 
random variables (resp. for probability measures). 

Theorem 1.1. Let {rij} C S be a subsequence of S satisfying Muj oo, such that the 
sequence {|)I))~(4)| : j > l} of non-negative numbers converges, that is: 

l/^n, (4)| ^ r], 

for some rj G [0,1]. Then 

In, A Mr,, (1.12) 

where Mr, was defined in ( 11 . 1111 . 

Since [KKW( IKW] showed that the^et of adherent points of {/ 2 „( 4 )}„g 5 is all of [—1,1], 
the result above clearly implies that Cn does not converge in distribution for Mn —>■ oo; 
in particular, if the sequence {|/lnj(4)|} does not converge, then the set of probability 
distributions associated with the random variables {Cnj} has at least two different adherent 
points in the topology of weak-* convergence. It would be desirable to formulate a uniform 
asymptotic result a la (| 1 . 12 p with no separation of the full sequence S into subsequences 
according to the angular distribution of A„ (still as Mn oo). This has two subtleties 
though. 

First, since there is no convergence in distribution, we need to couple the random variables 
on the same probability space and work with some metric on the space of probability 
measures; we choose to work with the L^-metrics, p G (0,2). Second, as, given a number 
n G 5, there is no limiting value rj of ^n(4:), for each n G S the candidate Mr, for the 
limiting random variable will bear 

V = Vn = |Mn(4)| 


rather than its limiting value. The following result is the desired refinement of Theorem 
11.11 Its proof (omitted) boils down to a standard adaptation of the proof of |Du( Theorem 
11.7.1], which is in turn an extension of the well-known Skorohod representation Theorem 
(see m Theorem 11.7.2]) to the framework of double sequences of probability measures. 

Theorem 1.2. On some auxiliary probability space (A,£/,¥) for every n G S there exists 
a coupling of the random variables Cn and A4|^^(4)| such that, as Mn —>■ oo, 




Cn-M 


\U'n ( 4)1 


PI 


0 , 


(1.13) 


for every p G (0,2), and ^ 

Tn ——)• 0, a.s. — P. (1-14) 

Relation (|1.14p is equivalent to saying that, for every sequence {nj} C S such that Mnj 
oo, P(/lnj — -^|/I„.( 4 )| —)• O) = 1. The fact that Theorem 11.21 is actually a strenghtening 
of Theorem 11.11 follows from the observation that, under the most natural coupling of the 
family of variables {A4^},^g[oq] we have 


E[|A4^i - A4r,2l] < c\r]i - ml, 
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for all r/i, 772 G [0,1] (with c > 0 an absolute constant). In fact, by the triangle inequality and 
an immediate computation, Theorem II.21 implies the stronger L^-convergence, p G (0, 2), to 
suitably coupled in ( 11 . 1211 . 


1.4 On the proofs of the main results 

In Proposition 13.21 we compute the Wiener-ltd chaos expansion for the nodal length Cn 
(m, i.e. a series converging in L^(P) of the form 

CO CO 

= ^WOi{^n\Cq) = '^Cn[q]- (1-15) 

g=0 q=0 


Here Cq, q = 0,1,... are the so-called ITiener chaoses (see M), namely the orthogonal 
components of the L^-space of those random variables that are functionals of some Gaussian 
white noise on T - while Cn{q] '■= proj(£n|C'g) denotes the orthogonal projection of Cn onto 
the g-th chaos. 

The decomposition (11.1511 is of independent interest, and entails in particular the vanishing 
of all the odd-order chaotic components and the term of order two, i.e. Cn\q\ = 0 if 
q = 2m -|- l,m = 0,1,... or q = 2. The precise analysis of the asymptotic behavior 
of the fourth-order projection in Proposition 12.21 will allow us to show that its variance is 
asymptotic to the total variance of the nodal length (see Proposition 12 . 3p : since the different 
components are orthogonal by construction, this will imply that all the projections other 
than the one on the fourth chaos are negligible. We notice that it is relatively easy to 
show that the contribution to the nodal length variance of each of the chaotic projections 
of order g 7 ^ 4 is negligible. It is in principle also possible to directly bound the total 
contribution to the variance of the sum of all these projections, thus establishing relation 
(11.71) independently. However, this task seems to be technically demanding, and would make 
our argument significantly longer. Since the asymptotic result (11.71) is already available from 
|KKW| . we do not pursue such a strategy in the present manuscript. 

As a consequence, to study the asymptotic behavior of Cn it will be sufficient to focus 
on the above-mentioned fourth-order component; Proposition 12.21 shows that along subse¬ 
quences {uj} satisfying the same hypothesis as in Theorem ll.il we have 


Cuj [4] 


^\av{CnM 


M 


r?) 


where Airi is as in ( 11 . 111 ) . 

We are then able to prove Theorem 11.11 thanks to Proposition 12.31 and Proposition 12.21 


1.5 Plan of the paper 

In ^2.1l we recall Wiener-Ito chaotic expansions, which we then exploit throughout the whole 
paper to prove the main results, given in ^ 2.21 ^ is devoted to the proof of the chaotic 
expansion for the nodal length fProposition 13.2)) . whereas in 21 we prove Proposition 12.21 
and Proposition 12.31 Finally, in 21 we collect the technical proofs of auxiliary lemmas for 
the results given in 21 
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2 Proofs of the main results 

The proofs of our results rely on a pervasive use of Wiener-ltd chaotic expansions for non¬ 
linear functionals of Gaussian fields; this notion is presented below in a form that is adapted 
to the random functions considered in the present paper (see e.g. [NPllPT] for an exhaustive 
discussion). 

2.1 Wiener Chaos 

Denote by {llfc}fc>o the usual Hermite polynomials on M. These are defined recursively as 
follows: Hq = 1, and, for k > 1, 


Hk{t)=tHk-i{t)-H'^_^it). 

Recall that HI := : k > 0} constitutes a complete orthonormal system in 

L‘^{R,m{R),-f{t)dt) :=L^( 7 ), 

where 7 (t) = (27r)“^/^e“*^/^ is the standard Gaussian density on the real line. 

The arithmetic random waves (jl.3jl considered in this work are a by-product of a family 
of complex-valued Gaussian random variables {a^ : A E defined on some probability 
space (D,.^,P) and satisfying the following properties: (a) every ax has the form xx + iy\, 
where xx and yx are two independent real-valued Gaussian random variables with mean 
zero and variance 1 / 2 ; (b) ax and a,- are stochastically independent whenever A ^ {r, —r}, 
and (c) ax = aZx- Define the space A to be the closure in L^(P) of all real finite linear 
combinations of random variables ^ of the form 

C = zax + za-x, 

where A E and z £ C. The space A is a real centered Gaussian Hilbert subspace of 
L2(p). 

Definition 2.1. For an integer q > 0 the g-th Wiener chaos associated with A, written 
Cq, is the closure in L^(P) of all real finite linear combinations of random variables of the 
form 

for A: > 1, where the integers pi,...,pk > 0 satisfy pi + ■ ■ ■ + Pk = Q, and is a 

standard real Gaussian vector extracted from A (note that, in particular, Cq = M). 
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Using the orthonormality and completeness of H in together with a standard 

monotone class argument (see e.g. |NPl Theorem 2.2.4]), it is not difficult to show that 
Cq _L Cm (where the orthogonality holds in the sense of L^(P)) for every q ^ m, and moreover 


L2(Q,a(A),P) = 0U,; 

q=0 

that is, every real-valued functional T of A can be (uniquely) represented in the form 

OO OO 

F = Y,Pm{F\Cq) = Y,F[q], ( 2 . 16 ) 

q=0 q=0 

where as before F[q] := proj(F | Cq) stands for the the projection onto Cq, and the series 
converges in L^(P). Plainly, F[0] = proj(F | Cq) = E[F]. 

A straightforward differentiation of the definition (11.31) of Tn yields, for j = 1,2 


2'J^i _^ 

djTnix) = ^ \jaxex{x), 


(2.17) 


(Ai,A2)eA„ 


(here dj = Hence the random fields Tn,diTn,d 2 Tn viewed as collections of Gaussian 

random variables indexed by x G T are all lying in A, i.e. for every x G T we have 

r„(x), diTn{x), d2Tn{x) G A. 


2.2 Proof of Theorem 11.11 

We apply the Wiener chaos decomposition (I2.16p on the nodal length 

OD 

Fn = ^ ^ 

q=0 


(2.18) 


in L^(P). The following proposition is a reformulation of Theorem o with the projection 
Cn [4] of the nodal length Cn onto the 4th order chaos replacing and it will be proven in 


Proposition 2.2. Let {nj} C S be a subsequence of S satisfying Muj oo, such that the 
sequence {|)I])~(4)| : j > l} of non-negative numbers converges, that is, 

ll^(4)| ^ V, 

for some r] G [0,1]. Then, the corresponding sequences of random variables converges in 
distribution to as defined in (jl.lip . that is. 


FnM 




Mr 


Moreover, 


Var(/:„,[4]) 


1 + Enj 

512 Ml 


(2.19) 


( 2 . 20 ) 














The next proposition, whose proof is given in ^4.21 entails that the fourth-order chaotic 
component gives the leading term in the expansion, i.e. its behaviour asymptotically dom¬ 
inates the nodal length on the torus. 

Proposition 2.3. For every {rij ■ j > 1} F S subsequence of S such that limj^oo A/’nj = oo 
and the sequence {|/I^(4)| : j > l} of non-negative numbers converges, 


Yax [Cn, - 


Equivalently, under the above assumptions we have that 

Var ~ Var (Tn^.[4]) . 


( 2 . 21 ) 


( 2 . 22 ) 


Proof of Theorem 11.11 assuming Proposition 12.21 and Proposition 12.31 The chaotic 
expansion ()2.18p and Proposition 12.31 entail that, as j -t-oo, 

[4] -|- Op(l), 

where op(l) denotes a sequence of random variables converging to zero in probability. Ac¬ 
tually, by linearity we have 

£n.[4] 

= (2.23) 


Y^Var(/:„^.) 


It hence follows that and the random variable Tn [4] have the same asymptotic dis¬ 
tribution. Proposition 12.31 together with (|2.19p and (12.2311 allow to conclude the proof, i.e. 
they immediately imply (jl.l2p . □ 

Remark 2.4 (On the length of u-level curves). For u E M, let us consider the total length of 
u-level curves for arithmetic random waves, i.e. the sequence of random variables {Cn-u}n£S 
defined as 

Cn-u := length(r“^{?x}). 

Of course, Cn-o = C-n- The behaviour of {Cn-u}n&s for u / 0 exhibits rather different 
characteristics than for the nodal case. Indeed, following a slightly modified version of 
the arguments we develop here (based on Green’s formula and the properties of Laplacian 
eigenfunctions - see i.e. [Roi §7.3 and p.l34]), it can be shown that the second order chaotic 
projection of is given by 


CnA‘A = 





TT , . . 2 1 




Mn 


Ed"'!’-!)- 


(2.24) 


AgAti 


where cf denotes the standard Gaussian density; note that (|2.24p confirms CnA = 0 in the 
nodal case (see ^1.41 and Proposition 13.21 (al). 

A few comments are in order. Let us first notice that the asymptotic variance of Cn-,uA\ 
satisfies, as n —)• -|-oo such that Mn —^ +oo, 


Yai{Cn-u[^]) 






8 AA 


(2.25) 
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The variance of the length of u-level curves for u ^ 0 can be derived exploiting the same 
computations as in [KKWj . and it is then possible to check the following asymptotic equiv¬ 
alence: for u / 0, as J\fn —^ + 00 , 


Var(£„;„) ~ Var(£„;„[2]). 

Hence, the variance of the length of non-zero level curves has a larger asymptotic order of 
magnitude than in the nodal case (compare ()2.25p - (j2.20p ): indeed, the former is dominated 
by the term corresponding to the second-order chaos, rather than the fourth. At u = 0, 
the second-order chaos component of the length of n-level curves vanishes exactly, and 
thus the variance has a lower asymptotic magnitude, consistently with the so-called Berry’s 
cancellation phenomenon [BellWl[IW2| . Also, because the second-order chaos term (|2.24p is 
proportional to a simple sum of independent, identically distributed, hnite-variance centred 
random variables (discounting repetitions coming from the symmetric structure of An), it 
is trivial to show that it exhibits limiting Gaussian behaviour, in marked contrast with the 
non-universal and non-Central Limit Theorem emerging in the nodal case. 

Remark 2.5 (On local statistics). Our method can be applied, in principle, to prove 
limit theorems for the nodal length within a proper subregion of the torus too. While the 
derivation of the L^-expansion into Hermite polynomials does not require any new ideas or 
techniques, some of the variance computations, and consequently the limiting distribution, 
will be affected. Note indeed that the main results of the present paper are obtained by 
exploiting some exact cancellations which are taking place when evaluating integrals of the 
eigenfunctions on the full torus. We hence leave these generalizations as a topic for future 
research. 


3 Chaotic expansion of jC^ 

In order to prove Proposition 12.21 and Proposition 12.31 we need to compute the Wiener- 
Ito chaotic expansion (j2.18jl of the random variable Cn] we refer to [KLj for analogous 
computations involving the length of level curves in the case of two-dimensional Gaussian 
fields on the Euclidean plane. 


3.1 Statement 

Let us introduce some more notation to properly state the main result of this section. The 
nodal length (jl.4p can be formally written as 


r — 

^n, — 


/ 5o{Tn{e))\\VTn{e)\\d9, 

JT 


(3.26) 


where do denotes the Dirac delta function and || 

Lemma 3.1] and ^3.2.ip . 

We shall often use the following easy result from [RW] : 
Lemma 3.1 ( |RW| . (4.1)). For j = 1,2 we have that 

YaildjTnix)] = ^ 
where the derivatives djTn{x) are as in srm . 


the Euclidean norm in (see [RW] 
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Accordingly, for x = (xi, X 2 ) G T and j = 1, 2, we will denote by djTn{x) the normalized 
derivative 

( 3 - 27 ) 

In view of convention (j3.27p . we formally rewrite (I3.26h as 


r — 

^T). - 


dvr^n 


(5o(Tli(x)) V dxTn{xY + d 2 Tn{xY dx. 


We also introduce two collections of coefficients {a 2 n, 2 m : n,m > 1} and {j32i : I > 0}, that 
are related to the (formal) Hermite expansions of the norm || • || in and the Dirac mass 
(lo(-) respectively. These are given by 

P 21 ■= —7^H2i{0), (3.28) 

V dvr 

where H 21 denotes the 2/-th Hermite polynomial, and 


^2n,2m 


1^ {2n)\{2m)\ 1 

V 2 n\m\ 2 n+mPn+m [jJ , 


where for N = 0,1, 2,... and x G M 


N 


Pn(x) := ■ (-1) 


N 


j=0 



(2j + 1)! 
0-!)^ 


(3.29) 


being the so-called swinging factorial restricted to odd indices. 

We are now ready to state the main result of this section. It illustrates the cancellations 
that occur for the components of the chaotic expansion (j2.18p of Cn (precisely, odd terms 
and the second-order one). Consistent to Proposition 12.31 computing the fourth-order com¬ 
ponent only is sufficient to establish the asymptotic behavior of the nodal length. However, 
we believe that the complete expansion is of clear independent interest; for instance, (a) it 
gives the basic building block to extend our results to other random fields on the torus and 
(b) it sheds some light on the Berry’s cancellation phenomenon [Bel IWll IW2| . as discussed 
also earlier in Eemark 12.41 

More precisely, as far as point (b) is concerned, we note that the nodal length Ci of 
Gaussian Laplace eigenfunctions Ti, £ G N, on the two-dimensional sphere have the same 
qualitative behavior. Indeed, on one hand in the chaotic expansion of Cn, the odd terms 
and the second chaotic projection vanish and the fourth-order component exhibits the same 
asymptotic variance as the full nodal length (see |R o| ). On the other hand, it is also shown 
in |Ro] that the second chaotic projection in the Wiener-Ito expansion of the length of 
level curves Tf^(u), n G M vanishes if and only if u = 0. These results explain why the 
asymptotic variance of the length of level curves is consistent to the natural scaling, except 
for the nodal case [Wl[ IW^ . Finally, we note that an analogous cancellation phenomenon 
occurs for the excursion area and the Euler-Poincare characteristic of excursion sets for 
spherical eigenfunctions, see |MW[ IMR] ICMW) . 

Proposition 3.2 (Chaotic expansion of Cn)- (a) For q = 2 or q = 2m-|-1 odd (m> 1), 


Cn[q] = 0, 
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that is, the corresponding chaotic projection vanishes. 
(b) For q >2 


^n[2g] 


47r^n 


q u 

EE 


0!-2k,2u-2khq-2u 

2 ^^^^{2k)\{2u-2k)\{2q-2u)\ 


(3.30) 


X / H 2 q- 2 u{Tn{x))H 2 k{dlTn{x))H 2 u- 2 k{d 2 Tn{x)) dx. 

Jy 

Consolidating the above, the Wiener-Ito chaotic expansion of Cn is 

Ot2k,2u-2kl32q-2u 


C-n — lEEn + 


. 9 +0O q u 

4Aeee 


/ 

Jl 


2 ^SfSfc;(2*=)K2..-2t)!(2,-2u)! 


X / H 2 q- 2 u{Tn{x))H 2 k{dlTn{x))H 2 u- 2 k{d 2 Tn{x))dx, 


in L^(P). 

3.2 Proof of Proposition 13.21 

Let us start with an approximating result in L^(P) for the nodal length 

3.2.1 Approximating the nodal length 

Consider the family of random variables > 0} defined as 

(3-31) 

where l[-e,e] is the indicator function of the interval [—e, e], and H-H is the standard Euclidean 
norm in 

In view of the convention (I3.27|l we rewrite (j3.31|) as 

M-e, 6 ]{Tn{x))^/difnixy + d2r„(x)2 dx. 

In [RWl Lemma 3.1] it was shown that, a.s. 

£„ = lim£^, (3.32) 

(a rigorous manifistation of (I3.26M . and moreover, by [RW( Lemma 3.2], is uniformly 
bounded, that is: 

< 12^/K. (3.33) 

Applying the Dominated Convergence Theorem to (I3.32|) while bearing in mind the uniform 
bound ()3.33l) implies that the convergence in (13.321) is in L^(P), i.e. the following result: 

Lemma 3.3. For every n £ S, we have 

iimE[|/:^-/:„|2] = 0. 


12 




























3.2.2 Proof of Proposition [3T^ technical computations 

In view of Lemma 13.31 we first compute the chaotic expansion of and then deduce 
Proposition 13.21 bv letting e ^ 0. Let us start by expanding the function into 

Hermite polynomials, as defined in fl2.ll 


Lemma 3.4. The following decomposition holds in (where, as before, 7 is the stan¬ 

dard Gaussian density on Mj.' 


1 1 

271[-m(') = E7TA'«'(■)> 

1=0 


Pi = -^7 (e) (Hi-i (e) - LI /-1 (-e)), 


/3o^ = 


2e 


7 (t) dt. 


where, for I > 1 

while for I = 0 

Moreover, as e ^ 0, 

Pf Pu 

where for odd £, Pi = 0 whereas Pi coincides with I13.28\} for even i. 

Proof. Using the completeness and orthonormality of the set H in one has that 


Po~'h f-£ ^ 

Now, if I is odd, then Hi-i is an even function, and therefore Pf = 0: it follows that 
1 1/1 \ 

= S + E 77 (-7^ (A A)j 


Using the notation (I3.28p . we have that, for all I > 0, 

limft', = -^(2i - = A,. (3.34) 

□ 

Note that the set {Pi : / = 0,1, 2,...} can be interpreted as the sequence of the coefficients 
appearing in the formal Hermite expansion of the Dirac mass 6q. 

Now hx X G T, and recall that the coordinates of the vector 


Vr„(x) := {diTn{x),d2T^{x)), 


are unit variance centered independent Gaussian random variables (see i.e., |KKW] L Now, 
since the random variable ||VT„(x)|| is square-integrable, it can be expanded into an (infi¬ 
nite) series of Hermite polynomials, as detailed in the following statement. 
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Lemma 3.5. For (^ 1 ,^ 2 ) a standard Gaussian bivariate vector, we have the -expansion 

+00 n 

11(2-^^)ll = E E 

n=0 m=0 ^ ' 

where the a 2 n, 2 n- 2 m as in (I3.29p . 

Proof. We may expand 

+00 U 


\\{z,,z2)\\ = y,Y 1 


^u,u—m 

^ u\(u — m)\ 
12=0 m=0 ' ' 


HuiZl)Hu-m{Z2), 


where 


0^77.n—m — 


hi ^ 


+ z^Hn{y)Hn-m{z)e 2 dydz. 


Our aim is to compute an,n-m as explicitly as possible. First of all, we observe that, if n 
or n — m is odd, then the above integral vanishes (since the two mappings 2 ; -s/y"^ -\- 

and y \/y‘^ F z^ are even). It follows therefore that 


H-cxd n 


\\iz,,z2)\ = j2J2 


^2n,2n—2m 

(2n)!(2n — 2m)! 

n =0 m =0 ^ ^ ^ ^ 


H2n{Zl)H2n-2m,{Z2). 


We are therefore left with the task of showing that the integrals 

— / \/y 2 + z‘^H2n{y)H2n-2m{z)e 2— dydz, 

JR2 


Q^2n,2n—2m — 


where n > 0 and m = 0,..., n, are given by (|3.29p . One elegant way for dealing with this 
task is to use the following Hermite polynomial expansion (see e.g. [SB Proposition 1.4.2]) 


+ 0O 






(3.35) 


a=0 


Let us consider the integral 


hi ^ 

4 vr Jr 2 


^ 2 ^26^2 dydz = 


if ^ 

4 vr Jr 2 


n (a-A )2 + (z -^)2 

^ + z^e 2 dydz. 


This integral coincides with the expected value of the random variable W := 
where (Y, Z) is a vector of independent Gaussian random variables with variance one and 
mean A and y, respectively. Note that W‘^ = has a non-central y^-distribution: 

more precisely, ~ X^(2, A^ + y?). Its density fiY 2 (see e.g. [ASl (26.4.25)]) is given 

by 

v2 I ,,2^ 




{t>0}) 


j=0 


where f 2 + 2 j is the density function of a X 2 -i- 2 j"bisti'ibuted random variable (see e.g. 
(26.4.1)]). Therefore, the density fw of IE is fw{t) = /iv 2 (t^) 2t, i.e. 


/„,(«) = g<(E±E)ZE/,^,^.(*2) 2i 1 


j=0 


Ht>o}- 
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Hence, the expected value of W is 

+ 00 


= r°° dt. (3.36) 

i=0 -^0 


From the definition of f2+2j in m (26.4.1)] we have 


r+cxD 


f2+2j{t )t dt = 


Substituting (|3.37l) into (13.361) we have 


r +00 

j) Jo 


2i+JT(l+i) 

21+jT(1+j) ■ 


t2J+2^-ty2dt 


(3.37) 


j=0 


E[1H] = g ((A^ +^^)/2y 


,fi). (3.38) 


Applying Newton’s binomial formula to ((A^ + fi‘^)/2y, we may expand the function F in 
p.38p as follows: 


F{X,fJ^) 2^ ^ ^ 

a=0 b=0 j=0 •’ 1=0 


E 


^ (-i)“A2“ ^ (-i)V“ 1 ^ , 2 ,2.-2ina(2i - i)x/| 


E ,1 E ( ; j 21+2ir{l + j) 


s ^ir C:!)^ s n.a(2i -1)^1 ^ /i\ 21+2 

9anl 3l92.7Tn m id Z^ I / / 


a 2j+2b-2Z 


2“a! 2^6! ^ 7'!22ir(l + 7) 

a,b=U i =0 ^ Z=0 ^ 2 

Setting n := ^ + a and m := j + b — I, we also have that 

FIA Ll) = ^ (-1)° (-1)^ ^ nz=l ~ i-)\/f f j\ ^2l+2a 2j+2b-2l 

^(A,iz) 2^ _ , 2^ 7 i 92 irn m 71 2^1/]'^ ^ 


2«a! 2^6! ^ 7'!22ir(l + j) 

a,b=0 7=0 \ ' JJ 


_y^.p n!=i(2i-l)vf ^ (-l)(r^-0 (_i)>»+z-7 /j 

Z^ Z^ j!22ir(l+j) 2’^“^(n —/)! 2™'+^“.?(m + ^ — j)! \Z 

Since F{X,/j) = E[W] from (|3.38p . on the other hand ()3.35ll yields 


A^’"^^”^.(3.39) 


F(A,/r) = 


1 


[ \/y2 + 2 "2 e ° 2 

i]R2 

1 7 \a ,.b 

iK. ^_n k_n 


/ 


+00 

E 

a,6=0 


-e 2 

\ 


1 


a!6!27r 7^2 


+ z‘^Ha{y)Hh{z)e t dydz 


V 


A“Ai^ (3.40) 


/ 


:=d(a,6) 

By the same reasoning as above, if a or 6 is odd, then d{a, b) in (I3.40p must vanish. By 
combining the expansions in (I3.39D and (l3.4Up . we have 


2 r _ 2+2^ 

a 2 n, 2 m = TT I V+ Z^H 2 n{y)H 2 m{z)e~^‘~^ dydz 
4 vr Jk2 

, (- 1 )-+- ^ (- 1 )^' n!=/( 2 i - 1 ) z’ 


= (2n)!(2m)!- 


2^2 H-m 


E 


E 


(i) 


2ij!r(l+j) ^ {n - l)\{m + I - j)\ 


-. (3.41) 
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The equality (j3.29l) now follows from (I3.4ip and some computations: 


^2n,2m — 


-/ 


+ z‘^H 2 n{y)H 2 m{z)e 2 dydz 


= (2n)!(2m)! 


(- 1 ) 


m+n 


nlU(2i - 


2 ^n-\-m 




w J 


E 


( 1 ) 


2ij!r(l+j) ^ {n - l)\{m + I - j)\ 


= (2n)l{2m] 


,(-l) 


m+n 


2n+m 


> 2j{j\y ^ 


(]) 


3 ' 1=0 

(2n)!(2m)! (-1++^ . (2j + l)!!+f ^ 

n!m! 2"+"^ ^ 2ij! ^ 


(n — /)!(m + / — j)! 


3=0 


1=0 


n 


m 




■(2n)!(2m) 

1 

n\m\ 

2n.+m 

■(2n)!(2m) 

! 

n\m\ 

2n.+m 

(2n)!(2m) 

! 

n\m\ 

2n,+m 


3=0 

n+m 

E(-i)^ 

j=0 

n+m 

i=o 


2Jj! 


(2(j + 1))! fn + m 


2^+^20j\{j + l)\\ j 


22i(j!)S 


□ 

We note that for two-dimensional random fields on the plane, the chaos decomposition 
of the length of level curves was derived earlier by Kratz and Leon, see m-, our derivation 
of the projection coefficients in Lemma 13.51 is different from theirs (albeit equivalent, by 
uniqueness), and it was hence reported for the sake of completeness. 

Proof of Proposition 13.21 In view of Definition 12.11 the computations in Lemma 13.41 
and Lemma [3. 5 1 (together with the fact that the three random variables Tn{x), diTn{x) and 
d 2 Tn{x) are stochastically independent, as recalled above) show that, for fixed a: G T, the 
projection of the random variable 


^l[-e,e\{Tn{x))\l diTn{xf + d2Tn{xf 

onto each odd chaos vanishes, whereas the projection onto the chaos C 2 q-, for g > 1, equals 

(J U p£ 

0!2m,2u—2mP2q—2u 


EE 


- 2m)!(2<, - 2ii)! 


H2q-2u{Tn{x))H2m{dlTn{x))H2u-2m{d2Tn{x)). 


Since fjdx < oo, standard arguments based on Jensen’s inequality and dominated conver¬ 
gence yield that C%[q\ = 0 if g is odd and for every g > 1, 


+ [2g] 




EE 


Ol2m,2u—2mP2q—2u 


2 ^^^^{2m)\{2u-2m)\{2q-2u)\ 


X / H 2 q- 2 u{Tn{x))H 2 m{dlTn{x))H 2 u- 2 m{d 2 Tn{x)) dx. 
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In view of Lemma [3^ and (12.181) one has that for every q > 0, as e —>• 0, necessarily 

converge to Cn[q] in We just proved that 


^n[<l] = 0 

for q = 2m + 1 as stated in part (a) of Proposition 13.21 Moreover, using (I3.34p . we deduce 
from this fact that representation (13.301) in part (b) of Proposition 13.21 is valid. To complete 
the proof of part (a) of Proposition 13.21 we need first to show that that Tn[2] = 0. From 
the previous discussion we deduce that T,i[2] equals 

f H 2 {d 2 Tn{x)) dx 
3 Jf 

Since H 2 {x) = — 1, we may write 

f H 2 {Tn{x))dx= f {Tn{xf-l) ^ OAdA'CA-A'(x) “ 1 j dx 

Jt Jt Jt J 

= T a,ayUy{x)dx-l = ±Y,iM'-^'>’ 

X^X'GAji y ^ y X^An 

sr 

where 5^' is the Kronecker symbol. (Observe that E[|aAp] = 1, hence the expected value of 
the integral H 2 {Tn{x)) dx is 0, as expected.) Analogously, for j = 1,2 we have 


(3.42) 


(3.43) 



H 2 {djTn{x)) dx = [ fV' XjX'jaxayex-\'{x) - l] dx 

hynNnx^x'eAx. J 




(3.44) 


where the used Lemma [Q to establish the last equality. 

Since a 2 n, 2 m = « 2 m, 2 n, and in light of (I3.43p and (|3.44l) . we may rewrite (|3.42l) as 


£„| 2 ] = 


= V dTT^ 


n / ao,o/52 1 /I 12 , “0,2/30 1 2 2 , \2',/| 12 

—w E (i“>i - 0 + —E 


Ac Ar? 


AgA„ 


V^^AT (^ “ 1 ) + 2q:o,2/3o ^ (|aAp - 1) 

y A^An AeAtx 

= VdTT^(“ 0 , 0^2 + 2ao,2/3o) ^ (I^Ap - !)■ 


Since 


1 Pk 


“0,0 — \ W , “0,2 — “2,0 — 


2 V 2 
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(3.45) 


























we have that 


«o,o/ 32 + 2ao,2/3o = 0 


and hence Cn\^] = 0 from (I3.45p . The proof of Proposition 13.21 is hence concluded, in view 
of (I2T8]) . □ 


4 Proofs of Proposition 12.21 and Proposition 12.3 


One of the main findings of the present paper is that, for any sequence {nj} such that 
Mrij —>■ oo and |/r„^(4)| converges, the distribution of the normalised sequence {Cn^} in 
(ll.lOp is asymptotic to one of its fourth-order chaotic projections. The aim of this section 
is a precise analysis of the asymptotic behavior of the sequence 


^rij [4] 

y/Var(£„^,[4])’ 


j > 1- 


which will allow us to prove Proposition 12.21 and Proposition 12.31 


4.1 Preliminary results 

Here we state the key tools for our proofs: first an explicit formula for [4] and then a 
Central Limit Theorem for some of its ingredients. First we need some intermediate results, 
whose proofs follow immediately from the fact that, for every n £ S, 

X] (Ai + A| - 6 A 1 A 2 ) , 

A=(Ai ,A2 )gA7i 


as well as from the identity Af -|-A| = n for A = (Ai, A 2 ) G A„, and the following elementary 
symmetry consideration: 


1 






1 


A=(Ai,A2)eA. 

Lemma 4.1. For every n £ S we have 

1 




^ ^2- 


nWn 




A—(Ai ,A2 )gA7i 


4 _ 3 -|- ^n(4) 


A—(Al ,X2)^A'n 


where £ = 1,2, and moreover 


1 


A=(Ai ,A2)£An 


2 \ 2 _ 1 A^n(4) 




Let us state now the above mentioned CLT result. Let us define, for n £ S 

( Wi{n) \ 


W{n) := 


W2{n) 

Wsin) 

\ w4n) y 


A2>0 


Z] (|oa|"-1) 


( n \ 

A? 


A2 

^2 


(4.46) 


\ A 1 A 2 / 

Exploiting the representation ()3.30l) in the case q = 2, one can show the following. 
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Lemma 4.2. We have, for diverging subsequenees {uj} C S sueh that Muj — ^ +oo and 
firiji^) converges, 

(l + - 2W2{n,f - 2W3in,)^ - 4W,{njf + op(l)) . (4.47) 

Uj 

The proof of Lemma 14.21 will be given in ^ 



Lemma 4.3. Assume that the subsequence {%•} C S is such that Mn^ +oo and /Inj(4) ^ 
r] E [—1,1]. Then, as nj —>■ oo, the following CLT holds: 


W{nj) 



( Z,\ 
Z2 
Zs 

\zj 


(4.48) 


where Z{rj) is a centered Gaussian vector with covariance 


S = S(,?) 


( 1 

1 
2 
1 
2 

V 0 


1 

3+?; 

A 

8 

0 


1 

iizz 

3+v 


0 

0 

0 

i-v 


0 V / 


(4.49) 


The eigenvalues of S are 0, and hence, in particular, S is singular. 

Proof. According to [NPt Theorem 6.2.3], in order to achieve the desired conclusion it is 
sufficient to prove the following relations: (a) for every fixed integer nj, each component 
of the vector W{nj) is an element of the second Wiener chaos associated with A (see 
112.11) . (b) as Uj —>■ oo, the covariance matrix of W{nj) converges to S, and (c) for every 
A: = 1, 2, 3,4, as nj —>■ oo, one has that Wk{nj) converges in distribution to a one-dimensional 
centered Gaussian random variable. Part (a) is trivially verified. Part (b) follows by a direct 
computation based on Lemma l4.11 as well as on the fact that the random variables in the 
set 

{joAp - 1 : A E Anj, As > 0} 

are centered, independent, identically distributed and with unit variance. To prove part 
(c), write A+ := {A E A„^., As > 0} and observe that, for every k and every Uj, the random 
variable Wk{nj) is of the form 


Wk{nj)= ^ Cfc(nj, A) X (IoaI^ - 1) 
AeA+ 


where {ck{nj,X)} is a collection of positive deterministic coefficients such that 


max Ck{nj, A) —)• 0, 
AeA+ 


as Uj ^ oo. An application of the Lindeberg criterion, e.g. in the quantitative form stated 
in [NPl Proposition 11.1.3], yields that Wk{nj) converges in distribution to a Gaussian 
random variable, and therefore that (c) also holds. This proves (I4.48p . Since it is easy to 
verify the claimed eigenvalues of S in (I4.49P via an explicit computation, this concludes the 
proof of Lemma 14.21 □ 
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4.2 Proof of Proposition [2T^ asymptotic behaviour of [4] 

Proof of Proposition 12.21 assuming Lemma I4.2L Let {rij} C S' be such that Afnj —)• oo 
and (4)| —>■ t/ G [0,1]. For each subsequence {n'j} C {rij} there exists a subsubsequence 
{n"} C {n' } such that it holds either (i) /i„»(4) rj or (ii) /r„//(4) ^ —r]. Set 


v(n^ := 




E„ 


512 AA2, 


j > 1 - 


Then, as n" 


oo 


Q(n') 


^n'J [4] 

u(n") 


-^1 + Zf- 2Zl - 2Zl - 4Z|, 


(4.50) 


by Lemma 14.31 and Lemma l4.2t here Z = ^( 77 ) G is as in (|4.48l) . i.e. a centred Gaussian 
4-variate vector with covariance matrix as in (I4.49p . 

Actually, the multidimensional CLT stated in (I4.48P implies that 


(VFi(n,f, W 2 {njf, Wsinjf, W^{njf) 


( ry^l ry2 ry^ v2'\ 


A simple computation of Gaussian moments now yields 
Var {Z\ - 2Zl - 2Z| - 4Z|) 


= 2 + ; 


3 + 77 


+ ' 


3 + 77 


+ 32 


1 — 77 


-2-2+4 


2,1 

= 77+1, 


entailing in particular that Var {Z\ — 2Z\ — 2Z| — 4Z|) is the same in both cases (i) - (ii). 
We can rewrite (I4.50p as, for tt," —)• + 00 , 


[4] 


1 


(1 + Zl - 2Zl - 2Z| - 4Z| 


+ 1+^ 77 ( 77 '') \/l + 772 

We claim that, in both cases (i)-(ii), the random variable 

-7=L=(i + Z2 - 2Z| - 2Z| - 4Z^ 

V 1 + ^ ^ 

has the same law as Air/, as defined in (|l.lip . 

To verify this, let S = 51 ( 77 ) be the covariance matrix of (Zi, Z 2 , Z 3 ), where Zj = 
for 7 = 1,2, 3, i.e. 


- 

VVarCZi) 


S(77) = 


We diagonalize S: 


/ 1 \/2 v'2 \ 

v'3+i? v'3+i? 

+2 1 1-v 

y/3-prj 3+?7 

]-2 1 

V V^+V 3+7? / 


S = ADA^ 
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where A is the orthogonal matrix 


A = 


0 


yv+ 7 

'3H^ 


\ 


/ ys+r] 

An+7 

v /2 _ _ 

Vv+T \/2 y2ArA^ 

\/2 _]_ y3+ri 

\ yrj+l \/2 \/ 2 \/» 7+7 / 


and 


Hence Zi = Xi, 


D = diag 


rj + l 2 ( 1 +r/) 
3 + r/’ 3 + 7? 


,0 


^2 = ^ 

y/^ + V V^ + V 

V _ w VI + ^ 

•^3 /q— j /q— i ^25 

+3 + 7? V3 + ?? 

where X = {Xi,X 2 ) G is a bivariate standard Gaussian random vector. Let us now set 

8 

Adding X^, as one more standard Gaussian, independent of {Xi,X 2 ), we find 


1 


\/i + V 


:(1 + zf - 2Zl - 2Zl - 4Z| 


+64V’2 - dSV’ + 10 


Xl - 27/; 




V3 + f? V3 + ?? 


- 2ip 


V^ Y _ V^ + ^ 

V3 + ?/ ^ V3 + ?? 


Xo 




1 


vT+V 

1 


-tt(3 + V 


- -(1 - 7?)X| + 1 

3 + 7? ^ 3 + 7? V 2^ ^ 


(2-(l + r?)X|-(l-r?)X|)^A4,. 

;Cn,[4] 


2 V 1 + ?/2 

In particular, we have therefore that 


K{nj) := 


■s/l + jf v{nj) 


Mr 


(4.51) 


To conclude the proof, we observe that, since {K[nj)'\ is a sequence of random variables 
belonging to a fixed Wiener chaos and converging in distribution, one has necessarily that 
fbv virtue e.e. of [NRl Lemma 2.1- (ii)]) 

supE|iL( 77 j)|^ < oo, Vp > 0. 

Uj 

Standard arguments based on uniform integrability yield therefore that, as rij —oo. 


Xai{Cr, 


v{njY{\ + r? 2 ) 
which is the same as (I2.20p . 


= nK{nj?) ^ nMl) = 1 , 


□ 
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4.3 Proof of Proposition [2.3^ C,nj[^ dominates 

Now we are able to prove one of the main findings in this paper, i.e. that the fourth-chaotic 
projection and the total nodal length have the same asymptotic behavior. 

Proof. Let us hrst prove ()2.22p . Note that (11.911 and Proposition 12.21 immediately give 
(j2.22p i.e., as —y -t-oo, 

Var(£„^.) ~ Var(£„^.[4]). 

Now, since different chaotic projections are orthogonal in L^, from part (b) of Proposition 
[Q we have 

+00 

YaiiCn,) = Var(/:„, [4]) + ^ Yai{Cn, (2q)). (4.52) 

q=3 

Dividing both sides of (14.521) by Var(£„^. [4]), we immediately conclude the proof of Propo¬ 
sition 12.31 

□ 


5 Proof of Lemma 14.2t explicit formula for [ 4 ] 

Consider the following representation of Cn^ [4], that is a particular case = 2 of (I3.30p : 

£„.[4] = J^Hi{Tn{x))dx (5.53) 

+ «o 4/3 o f j)^(a;)) dx + f H^idifnix)) dx + 

4! Jj 4! Jj 

+ X ^2(rn(x))/72(d2rn(x)) dx 

+ X ^2(rn(x))/72(dirn(x)) dx + 

+ X ^2(dirn(x))/72(d2r4x)) dx) , 

where the coefficients a.,, and fd. are defined according to equation (|3.29l) and equation 
(I3.28p . respectively. 

5.1 Auxiliary results 

The next four lemmas yield a useful representation for the six summands appearing on the 
right-hand side of (|5.53l) . In what follows, n G S and, moreover, to simplify the discussion 
we will sometimes use the shorthand 

E=E’E=E E=E’ 

A A=(Ai,A2)gA7t, A,A^ A,A^gA7t, A:A2>0 A=(Ai,A2)EAn 

A2>0 

in such a way that the exact value of the integer n will always be clear from the context. 
Also, the symbol {%•} will always denote a subsequence of integers contained in S such 

that Afuj —)• 00 and /inj( 4 ) —^ G [—1,1], as nj —>■ 00. As before, we write ‘—>■’ to denote 
convergence in probability, and we use the symbol op(l) to denote a sequence of random 
variables converging to zero in probability, as Mn —>• 00. 
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Following [KKW| . we will abundantly use the fine structure of the length-A spectral cor¬ 
relation set: 

‘S'nCd) ;= {(A, X', X", X'”) G : A + • • • + X'" = 0}. (5.54) 

Lemma 5.1 i |KKWj . p. 31). Let 5,^(4) he the length-A spectral correlation set defined in 
(j5.54p . Then 5,^(4) is the disjoint union 

Sn{^) = An{A)id Bn{A), 

where A„(4) is all the 3 permutations of 

I„(4) = {(A, A', -A, -A') : A, A' G A„, A / A'}, 

and Bn{A) is all the 3 permutations of 

Bn{A) = {{X,X,-X,-X) : Ag A,}. 

In particular, using the inclusion-exclusion principle, 

\Sn{4.)\=3J\fn{Afn-l)- 

Lemma 5.2. One has the following representation: 


f HfiTn{x))dx = 
JT 

Also, as Uj oo. 


6 




AA2 


^IuaI^ (5.55) 




6 . 


(5.56) 


O A 


Proof. Using the explicit expression H 4 {x) = x^ — 6x^ + 3, we deduce that 

[ Hi{Tn{x)) dx= [ (r„(x)^ - 6r„(x)2 + 3) dx 
JT Jt 

~ Hri cLfidyo-yo-x'" [ exp(27rz(A — A' + X” — X!”, x)) dx-\- 

-6-^ a^ay I exp(27ri(A — A', x)) dx + 3 


1 


A,A'eA„ 


jn Y1 axayayax'" . , 

^ A-A'+A''-A'''=0 ’■ 


1 


IoaI +3, 


(5.57) 


AgAti 


where the summation with the subscript A — A' + A" — A'" = 0 is over (A, — A', A", —X'") G 
5„(4). By the fine structure of SniAf described in Lemma IFiTl the right-hand side of (I5.57P 
simplifies to 

f HfiTn{x))dx = 3-^f ^ loAplaA'P - ^ IoaIA - 6-^ ^ loAp3 

"■ A,A'eA„ A "■ AeA„ 


AgAtj 


AgAti 


An \ x/Afn/‘^ 


^ (|«p-i) + <,Hi) -ttsE 


A:A2>0 


" A 


I^aI 
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where op(l) = 0 if is not an integer, otherwise 

0,(l) = (W,,/2rl/2(|„,„.„„,|2-l), 

thus yielding (I5.55P immediately. The limit (|5.56l) follows from a standard application of 
the law of large numbers to the sum, 


— y 

Mr ^ 


a\ = 


o A 


MnJ2 

"T A:A2>0 


|«a|^ + Op( 1 ), 


as well all the variables ax are i.i.d with 


E laAll =2. 


Lemma 5.3. For i = 1,2, 
24 


= E ^ (l“T -1)) + cp(i) -(f) 


M^ 


Mn y\A^A,A 2>0 

Moreover, as nj ^ oo, 

fij J J\rij ^ 

Proof. The proof is similar to that of Lemma 15.21 We have that 

[ Hi{difn{x)) dx = [ (difnix)^ - 6difn{x)^ + 3) dx 

Jr Jt 

= T ^—2 XeX'^X”X'”axdya\iiax"' / exp(27rz(A — A' + A" — A'", x)) dx + 

^ A,...,A'"eA„ 


□ 




1 4 
M^ rF 


A,A' 


Y 2 /* 

6—p —y XiX!faxd\i / exp(27ri(A — A', x)) dx + 3 

MnU^, Jj 

a\d\iaxiia\iii — ^ A^|aAp + 3 

AeAti, 

^ A^I«a|^ + 3 

\ A,A^ A / 


Ag Atj 


24 


f) 4E4 |«aI‘. (5.58) 


To conclude the proof, we first observe that the last term in the rhs of (|5.58p may be written 


as 




rij ) Mrij ^ 


= Op(l) + 




ni J Mn /2 

"T A:A 2>0 


(5.59) 


=:Ki{nj) 


=.K2{nj) 
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Now for the last term in the rhs of (I5.59P we have from Lemma l4.ll 


K2{nj) = 3(3 + /In(4)), 

so that the conclusion follows from the fact that /in(4) —>■ rj, as well as from the fact that, 
since the random variables {IoaI^—2 : A E Anj, X 2 > 0} in Ki{nj) arei.i.d., square-integrable 
and centered and < 1, EiLi(nj)^ = ^ 0. 

□ 


Lemma 5.4. One has that 

f H 2 {Tr,{x))(^H 2 {dlfn{x)) + H 2 {d 2 fn{x))') dx 


(5.60) 




AA2 

" A 


Proof. For ^ = 1,2, 


H 2 {Tn{x))H 2 {diTn{x)) dx = / {Tn{xf - l){diTn{xf - 1) dx 


I 1 7^yiaAaA'eA(a;)e_> 


)e-x'{x) “ 1 I I X] XiA”ax»ax'"ex»{x)e-y"{x) - 1 1 dx 


A", A'" 


2 1 


nJ\f^ 

" A-A'+A''-A'"=0 


E \ii \ III — — 

A(iA£ axax'ax'iax"! 


AL 


E 


ax - 


2 1 

nJ\f„ 


J^A||aAp + l. (5.61) 


An application of the inclusion-exclusion principle yields that the first summand in the rhs 
of (j5.6ip equals 


2 1 


nXf^ 

" A-A'+A''-A'"=0 


E \// \/// — — 

a^A£ a\a\ia\iiaxiii 


(5.62) 


nAL2 I 


+ 2 X! ^j^'j\ax\‘^\ax'\'^ — 




,A,A' A,A' 

Using the relation a_A = oa; we also infer that 


'^XjX'j\ax\‘^\ax'\‘^ = 


A,A' 


Substituting p5.63p into (I5.62p and then (I5.62p into p5.6ip we rewrite p5.6ip as 
[ H 2 iTn{x))H 2 {defn{x))dx 


(5.63) 




1 


nXf^ 

2 1 
nAfn 


A,A' 

E +1- 


Afn 


I«a| 


(5.64) 
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Summing the terms corresponding to di and 82 up, i.e. (I5.64p for £ = 1 and £ = 2, we 
deduce that the Ihs of (I5.60p equals 


H2{Tn{x))(^H2{diTr,{x)) + H2{d2Tn{x))') dx 


2 1 
nM^ 


(^1 + ^2) IoaI^IoA'I^ “ X] ("^1 ^2) 


.A,A' 


jr Y . E (^1 + ^ 2 ) l«Ap + 2 




n 


2 1 


nM^ 


, A,A' 


Mn 


nMri 


^ \A,A' A / "- A "a 

which equals to the rhs of (I5.60p . □ 

Our last lemma allows one to deal with the most challenging term appearing in p5.53p . 
Lemma 5.5. We have that 


H2{diTn)H2{d2Tn)dx 


= -4 


+4 


' i E Al(la.P-l) 


a/A/;/2 n 


A,A2>0 


- 4 


H 2 


E + 


A,A2>0 


n 2 


E (i«Ap-i)+op(i) 

V-^n/2 a,A2>0 


+16 


H 2 


1 1 


n 


E -^1-^2 (IoaI^ — 1) + op(l) 


^Mnl‘2^ A,A2>0 

and t/ie following eonvergence takes place as nj 00 : 


n2AA2 E^?^2 I«a| 


12 1 


n 


Proof. One has that 


2 A^2 
4 A 


4 1 


E^1^2|«a|^ a 3(1-7/). 


H2{diTn)H2{d2Tn)dx = E Ai A'l A'a'A'a'^d^aA' 


n2 A/"^ 

" A-A'+A''-A"'=0 

2 1 ,2 2 1 


OA" 


n A/), 


E^+aP 


n A/), 


E ^ikAp + 1- 


(5.65) 


(5.66) 


26 





























First of all, we note that for the first two terms in (I5.66P 


E 




= E 


A/" 


= 2. 


Let us now focus on (j5.66p . Using the structure of 54(n) in Lemma ETJ we obtain 

" A-A'+A''-A"'=0 


(5.67) 


4 1 


n- 


E'^i^'^2)^I«aP|oa'P + 2^ AiA2A^A2|aAp|aA'P - 3^ AfA||aA|^ 

A,A' A,A' A 


Let us now denote 


A := 


4 1 


n- 


EA?(A'2)>a|>A'P, 


A,A' 


^ ■“ AiA2AiA2|aAp|aA'P, 

^ A,A' 

A 


D := 




SO that (|5.66l) with (15.671) read 


(5.68) 

(5.69) 


H2{diTn)H2{d2Tn) dx = A + B + C + D. 


(5.70) 


We have that 


4 1 

I 

A,A' 

4 1 1 

4 1 1 

^A^2 


E^l('^2)^l«A|^|aA'|^ + E-^i(-^2)^|Qa|^|QA'|' 


A,A' 


A,A' 


E(« - Ai)(A'2)2|aA|EA'P + E “ (^i)^)|aAp|aA'P \ , 

A,A' A,A' 
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which we may rewrite as 


4 1 1 




A,A' 


A,A' 


I ^IZ(^2)^|aA|VA'P + n^A?|aA|>A'P 

"■ [ A,A' A,A' 

I - E Ai(A'2f |„j|2|a,,|2 - ^ A;(A;)2|ai| V,,|= 

A, 


A,A' 


n‘^X^2 

From ()5.69D and (15.711) . we get 
A 4 11 

^ + ^ “ ~^m2 


{ ^X](^2)^|aAlVA'P + n^A?|aAp|aA/|^ 

A,A' A,A' 


(5.71) 


1 2 


4 1 1 


+ 


4 1 
^2 


= -4 


-4 


^A^doAp - 1) 

L A 

E(I“aI"-i) 

A 


1 2 


^AfdttAp - 1) 


L A 


n 2 


(5.72) 


n 2 


1 1 


VKj 2 


n 


Y. A?(|o,p-l)+o,(l) 


A,A2>0 


+4 


H 2 


-7^ E (l»Ar^-l) + or{l) 

y-Nn/2 ^ ;i^2>0 


On the other hand, 


4 1 


4 1 


^ „2^22Z]^i^2A'iA2|aAp|aA'P ^2^2^ 

A,A^ 


n 2 


y~! A 1 A 2 (ioai^ - 1) 


= 16 


L A 

n 2 


1 1 




n 


y~l A1A2 (IoaI^ - 1) + op(l) 


A,A2>0 


(5.73) 


The first statement of Lemma 15.51 then follows upon substituting (15.721) , (15.731) and (15.681) 
into (15.701) . 
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Now to prove ()5.65l) it suffices to write 


12 1 


Y,XlXl\a,\^ 


^ - 2) + —2 W X] 






(5.74) 


-^^i^) + ~.Trr ^ ^i^2(I«a|^ - 2) + ;^-^^AiA2, 


rii^ Af^. 


^ \:\2>0 


VAA2. ^ 


and then Lemma [4.1l and an argument similar to the one that concluded the proof of Lemma 
15.31 allow to prove the result. 

□ 


5.2 Proof of Lemma 14.21 technical computations 

Substituting the results of Lemmas 15.2115.51 into (15.531) we obtain, as nj +oo, 



(5.75) 


where R{nj) is a sequence of random variables converging in probability to some constant 
G M. Computing the coefficients a.^. (see (|3.291) 1 and (5. (see (|3.281) 1. from (I5.75P we obtain 
that, as rij +oo, 

^n, [4] = J (vLi(n,)2 - 2 W 2 {n,f - 2W^{n,f - + i?(n,)) , 

where hLfc(nj), k = 1,2,3,4 are as in (|4.46p and R{nj) is a sequence of random variables 
converging in probability to 1. The proof is now complete. 

□ 
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